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Abstract 

The Holstein Hubbard and Holstein t-J models are studied for a wide 
range of phonon frequencies, electron-electron and electron-phonon inter- 
action strengths on finite lattices with up to ten sites by means of direct 
Lanczos diagonalization. Previously the necessary truncation of the phononic 
Hilbert space caused serious limitations to either very small systems (four 
or even two sites) or to weak electron-phonon coupling, in particular in the 
adiabatic regime. Using parallel computers we were able to investigate the 
transition from 'large' to 'small' polarons in detail. By resolving the low-lying 
eigenstates of the Hamiltonian and by calculating the spectral function we can 
identify a polaron band in the strong-coupling case, whose dispersion deviates 
from the free-particle dispersion at low and intermediate phonon frequencies. 
For two electrons (holes) we establish the existence of bipolaronic states and 
discuss the formation of a bipolaron band. For the 2D Holstein t-J model 
we demonstrate that the formation of hole-polarons is favoured by strong 
Coulomb correlations. Analyzing the hole-hole correlation functions we find 
that hole binding is enhanced as a dynamical effect of the electron-phonon 
interaction. 



PACS number(s):: 71.27.+a, 71.38.+i, 74.25.Kc, 75.10.Lp 
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I. INTRODUCTION 



Following the discovery of high-temperature superconductivity in the ceramic copper 
oxides, novel purely electronic pairing mechanisms due to the strong Coulomb correlations 
within the Cu02 planes have been investigated in detail. Recently, however, it has become 
clear that the lattice degrees of freedom are essential in understanding the puzzling normal- 
state properties of the cuprates [|T]-0]. Even if it should turn out that the electron-phonon 
(EP) interaction is not the relevant pairing interaction in those materials, its effects need to 
be reconsidered for the case of strong electron-electron interactions and low effective dimen- 
sionality as realized in the high-Tc superconductors. In particular, polaronic effects are sug- 
gested to play a non-negligible role in the copper-based materials La2-xSrxCu04+y [^-^ and 
even more in the isostructural nickel-based charge-transfer oxides La2_xSrxNi04+y [|7|,p!0 



Experimentally, photo- induced absorption experiments infrared spectroscopy [0 as 



well as infrared reflectivity measurements [jT3| unambiguously indicate the formation of 'self- 
localized' polaronic states (small polarons) in the insulating parent compounds La2Cu04+y 
and Nd2Cu04_y of the hole- and electron-doped superconductors La2_xSrxCu04+y and 
Nd2_xCexCu04_y, respectively. Therefore a growing theoretical interest in the study of 



strongly correlated EP models can be found in the recent literature [p!^-p6|. 

Probably the simplest microscopic models including both the electron and phonon de- 
grees of freedom are the Holstein Hubbard model 

Hh-h = -tYl {4aCja + H.c.) +UJ2 ^i^^il 

{ij)a i 

-y^^(6l + 6Jn, + huJ2{blb, + l) (1) 
i i ^ 

and the Holstein t-J model 

HH-t-j = -tY. (claCjv + H.c.) + [SiSj - ^n^nj') 

(bl + h) K + huY. (bh + 1) , (2) 
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where cjP annihilates (creates) an electron at Wannier site i with spin projection a, 
ni = rii-^ + riii, ^^"^ ^ denotes the transfer amplitude between nearest-neighbour (NN) 
pairs (ij). 'Hu-t-j acts in a projected Hilbert space without double occupancy, i.e., 
cj? = - n^^), and Si = \ Eaa' ct^aa'^ia' ■ The first two terms in (|I]) and (g) represent 
the standard Hubbard model and t-J model, respectively, where U is the on-site Coulomb 
repulsion and J measures the NN antiferromagnetic exchange interaction strength. The third 
and fourth terms take into account the EP interaction and the phonon energy in a harmonic 
approximation. Here, the on-site electron (hole) occupation number Ui {hi = 1 — hi) is 
locally coupled to a dispersionsless optical phonon mode, where Sp is the EP coupling con- 
stant, uj denotes the bare phonon frequency, and 6^ are the phonon annihilation (creation) 
operators. In the context of an effective single-band description of the copper/nickel oxides. 



the collective Holstein-coordinates = ^hjlMuj (&| + 6j) may be thought of as representing 
an internal optical degree of freedom of the lattice site i.e., in this case the dominant source 
of EP coupling is assumed to result from the interaction of dopant-induced charge carriers 
with the apical out-of plane or the bond-parallel in-plane breathing-type displacements of 
oxygen atoms. 

Unfortunately, for strongly coupled EP systems exact results exist only in a few spe- 
cial cases and limits p7|^0|. Whereas, in an approximative treatment, the weak-coupling 
regime [SpIt <^ 1) is well understood and dealt with by perturbation theory, the standard 
strong-coupling Migdal-Eliashberg theory |]3ll , |32| based on the adiabatic Migdal theorem 
might break down for strong enough EP interactions [Spji ^ 1) due to the familiar pola- 
ronic band collapse [^. Note that in the presence of strong Coulomb correlations, a rather 
moderate EP can cause a substantial reduction of the coherent band motion making the 
particles susceptible to 'self-trapping' |jT6|,^. The (single) polaron problem has been tack- 
led in the strong-coupling adiabatic {Tiuj/t <ti 1) and antiadiabatic {Tiuj/t ^ 1) limits using 
the Holstein |Q and Lang-Firsov 0] approximations, respectively. Both approaches yield 
a narrow polaronic band with an exponentially reduced half-bandwidth ||18[. Whether these 
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small polarons (or bipolarons) can exist as itinerant band states is still a heavily debated 



issue ||T|. Apart from variational calculations p5|- |37| , p0| little is known for intermediate 
values of EP coupling and phonon frequency Ep ^ huj ^ t and, in particular, for the many- 
polaron problem. In principle, exact diagonalization (ED) p^j38|JT8|j26[| and (quantum) 
Monte Carlo pPj-^Hf methods including the full quantum nature of phonons can close this 
gap. However, by using direct ED techniques it is necessary to truncate the phononic Hilbert 
space, and hence the accessible parameter space is limited by the size of the matrix one can 
diagonalize. Therefore ED studies up to now were limited to either small values of Ep, to the 
so-called frozen phonon approximation p4|-|47|], or to very small systems . In a 

previous work p4|, the authors have proposed a variational Lanczos diagonalization tech- 
nique on the basis of an inhomogeneous modified variational Lang-Firsov transformation 
(IMVLF) that allows for the description of static displacement field, polaron and squeezing 
effects in terms of the Holstein t-J and Holstein Hubbard models on fairly large clusters. 
Although the adiabatic and antiadiabatic as well as the weak- and strong-coupling limiting 
cases are well reproduced in this approach, the situation becomes less favourable at inter- 
mediate EP couplings and phonon frequencies and, in particular, in the crossover region 
from large-size nearly free polarons (FP) to small-size 'quasi-localized' polarons (i.e., in the 
vicinity of the so-called 'self-trapping' transition). Obviously, this regime requires a more 
accurate treatment of the phonons as quantum mechanical objects. 

Encouraged by this situation it is the aim of the present paper to perform a direct 
Lanczos diagonalization of the Holstein Hubbard and Holstein t-J models, preserving the 
full dynamics of quantum phonons. In particular, we investigate for the first time the 
low-lying excitations (spectral functions) on large enough lattices, in order to identify the 
dispersion relation of the (bi)polaronic quasiparticles. 
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II. COMPUTATIONAL PROCEDURE 



A general state of the model Hamiltonian Hh-h [Hh-i-j] describing Nei = + Ni 
electrons on a finite D-dimensional hypercubic lattice with N sites can be written as the 
direct product 

m=T.^f\l)el®\k)pk, (3) 
l,k 

where / and k label the basic states of the electronic and phononic Hilbert space with 
dimensions D^i = (^^^^ D^i = (^^^^ (^jvf^^) ^"^^ ^p^^ ~ ^respectively. Since the 

bosonic part of the Hilbert space is infinite dimensional we use a truncation procedure |48 
restricting ourselves to phononic states with at most M phonons: 

i^)p'^=n^(^i)"' io)p/. (4) 

i=i Jn'-l 

with 

j:n1<M, (5) 

i=l 

and 1 < A; < D^^^ = (M + A^)!/M!A^! . To further reduce the dimension of the Hilbert 
space, in the case of TLh-h we separate out the center of mass motion by transforming 
to new phonon operators Bf\ which can be taken into account analytically as displaced 
harmonic oscillators. For the Holstein t-J model it is more effective to exploit the point 
group symmetries of the original basis (^. 

Then the resulting Hamiltonian matrix is diagonalized using a standard Lanczos 
method. As the convergence of the Lanczos procedure depends on the (relative) differ- 
ence of neighbouring eigenvalues, l-Ej+i — Ei\/\Ei\, one needs to be very careful in resolving 
eigenvalues within the extremely narrow small-polaron band. To monitor the convergence 
of our truncation procedure as a function of M we calculate the weight of the m-phonon 
states in the ground state jip'o) of 

AT 









:icfi 



where m = ^ . (6) 

,k 1=1 
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At fixed M, the curve |c'"p(m) is bell-shaped and its maximum corresponds to the most 
probable number of phonon quanta in the ground state. To illustrate the M-dependences 
of the ground-state energy and the coefficients jc™"!^, we have shown both quantities for 
the single-electron Holstein model in Fig. |1|. In the numerical work convergence is achieved 
if the relative error of the ground-state energy is less than 10"''. In addition, we check that 
Eq is smaller than the estimate obtained from the IMVLF-Lanczos treatment of the phonon 



subsystem 



We have written the program in Fortran90 and ran it on a 64-node CMS. We were able 
to diagonalize Hamiltonian matrices up to a total dimension {Dtot) of about 82 millions. 
Since a matrix vector multiplication for this matrix size takes less than 150 seconds, the 
limiting factor of our numerical algorithm is the available storage. 



III. NUMERICAL RESULTS 



A. Holstein Hubbard model 



1. One-electron case 



In the first place, we investigate the polaron properties of the Holstein model with a 
single electron on finite lattices with up to ten sites using periodic boundary conditions. 
In the light of the literature over at least the last two decades p9| , |39| , p6| , p8| -pO[1 we expect 
a gradual transition from a (nearly free) large-polaron solution to a small-polaron-like 
ground state upon increasing the EP coupling. Since, in particular in the adiabatic regime, 
the formation of a polaronic state is accompanied by a strong reduction of the coherent 
electron motion, this effect should be observable in the expectation value of the kinetic 
energy Ep^kin/t = -E<ij>^{%\{ '^la'^ja ~^ H.c.) 1*^0)5 where \^o) is the ground— state wave- 
function. We therefore define an effective polaronic transfer amplitude p4| , 



U)/Ep^k^n{0,U), (7) 
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in order to characterize the increase in the quasiparticle mass |2J]. Note that tp,e// sub- 
stantially differs from the (exponential) polaron band renormalization factor (p) obtained 
analytically in the non-adiabatic Lang-Firsov and adiabatic Holstein cases . 



We illustrate the dependence of this effective hopping amplitude on the EP interaction 
strength in Fig. ^ where we have plotted tp,e// as a function of Ep at different phonon 
frequencies. First it is important to realize that there are two complementary (adiabatic 
and non-adiabatic) regimes for the polaronic motion. In the non-adiabatic regime, where 
the lattice fluctuations are fast and the phonons are able to follow immediately the electronic 
motion forming a non-adiabatic Lang-Firsov polaron (NLFP), one observes a very gradual 
decrease of tp,e// as Ep increases. At the same time the 'phonon distribution function', Ic"^]^, 
gets wider but the maximum is still located at the zero-phonon state. In the adiabatic 
regime, one notices a crossover from a large-size polaron (LP) in ID or nearly free polaron 
(FP) in 2D, described by a tp,e// that is only weakly reduced from its noninteracting value, 
to a less mobile (small-size) adiabatic Holstein polaron (AHP) for large Ep. We point out 
that the nature of 'delocalized' polaronic states, occurring in the weak-coupling region, 
is different in ID and 2D In the ID case, the FP state becomes unstable at any 

finite EP coupling. As expected the transition to the AHP state occurs if the EP coupling 
approximately exceeds half the bare electronic bandwidth and, in accordance with Monte 



Carlo results [p9| , ^ , is much sharper in two dimensions (in the remainder of this section 
we focus on the ID case). Nonetheless, all physical quantities are smooth functions of Ep, 
in particular there are no ground-state level crossings, i.e., the transition from LP/FP to 
AHP is continuous and not accompanied by any non-analyticities. While in the weak- 
coupling case we have rrimax = and the inclusion of higher phonon states (m > 5) does 
not improve the ground-state energy at all, in the adiabatic strong-coupling case {Ep = 4, 
huj = 0.4), the maximum in |c™p is shifted to multi-phonon states {rrimax — 8) and we need 
about 16 phonons to reach a sufficient accuracy within our truncation procedure. Note that 
a similar behaviour can be observed in the non-adiabatic regime {hu > t) provided that 
Ep ^ hu, e.g., for hu = 3 and Ep = 8 {Ep = 10) we find nimax — 2 in ID (2D). These results 
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confirm previous findings for tlie Holstein Hubbard model on very small size clusters (with 
two or three sites), where, as Ep increases in the adiabatic regime, a strong increase of the 
average number of phonons, (Nph), contained in the ground state, was observed (cf. Tab. I 
in Ref. Jl^ and Tab. I in Ref. |15|). In the center of mass system, the phonon expectation 



value in the polaronic ground state may be derived from the phonon distribution function 

To elucidate the difference between the 'extended' LP and 'quasi-localized' AHP states 
in more detail, we have calculated the electron-phonon density correlation function 

C,l.p,{\^ ~ j\) = {%Kb]b^\^o) , (8) 

which measures the correlation between the electron occupying site i and the density of 
phonons on site j [^. Results for Cei^ph{\i — plotted in Fig. |] at /ic<j = 0.4 for all 



distances i — j := Ri — Rj , show that for small Ep the correlation between the electron and 
the phonons is pretty weak and exhibits little structure, i.e., the few phonons contained in 
the ground state are nearly uniformly distributed over the whole lattice. In contrast, in the 
case of large EP coupling {Ep = 3), the phonons are strongly correlated with the position 
of the electron, thus implying a very small radius of the polaron. Note, however, that the 
translational invariance of the ground state is not broken. Since a polaron's mass is inversely 
proportional to its size, the AHP formed at large Ep is an extremely heavy quasiparticle. 
As can be seen from the inset of Fig. ^ the on-site electron-phonon correlation increases 
dramatically around the same value of Ep at which tp,e// becomes depressed (cf. Fig. ^. This 
means, in the adiabatic regime a strong short-range EP interaction can lower the energy of 
the system due to a deformation-potential-like contribution sufficiently to overcompensate 
the loss of kinetic energy. Nonetheless, the 'quasi-localized' (self-trapped) polaronic state 
has band-like character, i.e., the AHP can move itinerantly. 

In order to discuss the formation of a small-polaron band one has to calculate the 
low-lying excited states. As a first step, in Fig. ^ we classify the lowest eigenvalues of the 
Holstein model according to the allowed wave-vectors of the eight-site lattice for various 



phonon frequencies at Ep = 3. Here the 'band dispersion' Ek — Eq is scaled with respect to 
the so-called coherent bandwidth AE = sup^- E^ — infi^ E^- AE strongly depends on both 
ratios Sp/Tiuj and Ep/t, for example, we found AE{ep = 3, hu) = 0.0157, 0.1957, 2.9165, and 
4.0 for huj = 0.4, 0.8, 10.0 and oo, respectively. Of course, the simple Lang-Firsov formula, 
AElp = AD exp[—ep/huj], gives a good estimate of the polaronic bandwidth only in the non- 
adiabatic regime: AEipi^p = 3, hu) = 0.0022 {hu = 0.4), 0.0941 (0.8), 2.9633 (10.0), 4 (oo). 
Besides the strong renormalization of the bandwidth in the low-frequency strong-coupling 
regime it is interesting to note that the deviation of the polaron band dispersion from a 
(rescaled) 'cosine-dispersion' of noninteracting electrons is most pronounced at intermediate 
phonon frequencies hu ~ t, i.e., in between the extreme adiabatic (AHP) and antiadiabatic 
(NLFP) limits. This deviation may be due to a residual polaron-phonon interaction, with 
the phonons sitting on sites other than the polaron. To demonstrate that the low-lying 
eigenvalues do indeed form a well-separated quasiparticle band in the adiabatic strong- 
coupling regime {Sp = 3, hu = 0.4), in the inset of Fig. ^ we have displayed the lowest few 
eigenvalues in dependence on Ep. In the very weak-coupling regime {Sp = 0.5) the eigenvalues 
are barely changed from their Ep = values, where additional eigenvalues, separated from 
the ground-state energy Eq by multiples of hu (e.g., E2, -E3, and E4), enter the spectrum. As 
Ep increases a band of states separates from the rest of the spectrum. These states become 
very close in energy and a narrow well-separated energy band evolves in the strong-coupling 
case {Ep = 3). Obviously, the gap to the next higher band of eigenvalues is of the order of 
the bare phonon frequency hu. Neglecting degeneracies one may tentatively identify those 
five states as the states of the small-polaron band on the eight-site lattice. 

Keeping this identification in mind, in Fig. ^ we have plotted the lowest eigenvalues as 
a function of the (ID) i^-vectors belonging to various system sizes (A^ = 6, 8, 10). One 
notices that the dispersion E^ is rather size independent, i.e., the E^ values obtained for 
larger systems just fill the gaps. Undoubtedly, the smooth shape of Ek already provides 
good reasons for a quasiparticle band description of the AHP in the strong-coupling regime. 
To further substantiate this quasiparticle interpretation, we also have calculated the one- 
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particle spectral functions 

n 

with Nei = 1 for the non-equivalent i^-values of the six-site system using a polynomial 
moment method |^ . The idea is to see a direct verification of the coherent band dispersion 



Ex in terms of Ax{E). The electronic spectral functions Ax{E) are shown in the four insets 
of Fig. 1^. The important point, we would like to emphasize, is that the position of the first 
peak in each spectral function Ak{E) exactly coincides with the corresponding i?i^-value 
and the other peaks are at higher energies than any of the coherent band-energy values. 
This means, our exact results for the low-energy excitation spectrum of a single electron 
corroborate the existence of heavily dressed polaronic quasiparticles, where the electronic 
and phononic degrees of freedom are strongly mixed. Of course, in the very high-energy 
regime the results for Ak{E) can not be trusted just due to the errors induced by the 
necessary truncation of the phononic Hilbert space. 



2. Two-electron case 

Next, we wish to discuss the two-electron problem. Here it is of special interest to 
understand in detail the conditions under which the two electrons form a bipolaron. Whether 
or not a transition to a bipolaronic state will occur depends sensitively on the competition 
between the short-ranged phonon-mediated, i.e., retarded {hu < oo), attraction (oc Sp) and 
the instantaneous on-site Hubbard repulsion (oc U). 

We start again with a discussion of the mobility of the particles. Fig. ^ (a) shows the 
strong (gradual) reduction of the effective polaronic transfer amplitude tp,e// as Sp increases 
in the adiabatic (non-adiabatic) regime. Now let us mainly focus on the physically more 
interesting regime of 'small' phonon frequencies, hu = 0.4. In the case of vanishing Coulomb 
interaction U = any finite EP interaction causes an effective on-site attraction between 
the electrons forming a bipolaronic bound state (remember, e.g., that t/g// = U — 2ep follows 
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from the simple Lang-Firsov approach). This means, in the pure Holstein model the state 



with two nearly free (large) polarons does not exist, at least in one spatial dimension 
In the weak-coupling limit, the two-polaron state can, however, be stabilized by taking into 
account the on-site Coulomb repulsion. In this case, a crossover from a state of two mobile 
large polarons to an extended bipolaronic state occurs. The 'transition' will be shifted to 
larger EP couplings as U increases (see Fig. |] (a)). For example, aX U = 6 and hu = 0.4 
(3.0), we find that the binding energy of two electrons, E"^ = Eo{2) — 2Eq{1), becomes 
negative at about Sp = 1.7 (2.8). Further justification for this interpretation can be found 
from the behaviour of the effective bipolaronic transfer amplitude ||53|| , 



tb,eff = Eh^kiniep, t/)/-Ef,^fcm(0, U) (lO) 

with 

EbM^p, U)/t = - ^ {%{ep, U)\{c\^clc^^c^^ + H.c.)|iZ>o(£p, U)) , (11) 

<ii> 

shown in Fig. ^ (b). tfe.e// describes the coherent hopping of a on-site bipolaron from site i to 
site j. Contrary to tp,e//, at low EP coupling strengths, the bipolaronic hopping amplitude 
^6,e// grows with increasing Sp showing the increasing importance of the correlated motion 
of two electrons (but, quite clearly, we have \Eh^kin\ < \Ep^kin\)- At large EP couplings (e.g., 
for ep>l at U = and huj = 0.4), the on-site bipolaron becomes more and more localized 
and accordingly we observe a drop in tb^eff which corresponds to the drop in tp.e// in the 
case of one electron at the parameter values where the AHP becomes stable. Hence we will 
call this quasiparticle an adiabatic Holstein bipolaron (AHBP). 

To better illustrate the effect of pair formation in the ID Holstein (Hubbard) model, 
we present in Fig. ^ the electron-electron density correlation function 

Cei-eii\^ ~ J\) = {%iep,U)\n,n^\%iep,U)) - {%iO, U)\n,n^\UO,U)) (12) 

in the adiabatic regime with (b) and without (a) Hubbard repulsion. In each case we have 
displayed the results for C^i_^i{\i — j|) as a function of Sp in comparison to the electron- 
phonon correlation function Cg,_p^(|i — j\) given by (||). As Fig. |^ (a) shows, in the limit of 
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vanishing Coulomb interaction the on-site electron-electron correlation Cei-ei{0) dominates 
the inter-site correlations Cei-ei{\i — j\) with \i — j\ > 1, in particular for ep>0.9, i.e., 
in the AHBP regime where both electrons are mainly confined to the same site sharing 
a common lattice distortion. Therefore the transition from a mobile large bipolaron to a 
'quasi-self-trapped' on-site AHBP is manifest in a strongly enhanced Cei-ph{0) as well (see 
inset). Moreover, the transition should be associated with a significant reduction of the local 
magnetic moment, mioc{£p, U) oc (ip'oK^iT ~ '"^■ilYWo)) indicating the local pairing of spin up 
and down electrons. Indeed we found mioc{£p = l)/mioc{£p = 0.9)|^=o,7ia;=o.4 = 0.66. As can 
be seen from Fig. ^ (b), a somewhat different scenario emerges in the presence of a finite 
Coulomb interaction. Here, the Hubbard repulsion prevents the formation of an on-site 
bipolaronic bound state in the weak EP coupling regime. On the other hand, as recently 



pointed out by Marsiglio p6[, the retardation effect of the EP interaction may favour the 
formation of more extended pairs. That is, due to the time-delay the second electron can 
take the advantage of the lattice distortion left by the first one still avoiding the direct 
Coulomb repulsion. In fact, increasing the EP interaction, we find that both the nearest- 
neighbour electron-electron and electron-phonon density correlations starts to rise, while 
the on-site correlations remain small (cf. Fig. |^ (b)). Consequently, we may label this state 
an adiabatic inter-site bipolaron. We expect that at larger values of Ep the short-range EP 
interaction overcomes the Hubbard repulsion and as a result the two electrons coalesce on a 
single site forming a 'self-trapped' bipolaron. Unfortunately we are unable to increase the 
dimension of the Hilbert space to contain a large enough number of phonons in the adiabatic 
very strong-coupling regime. 

As already mentioned for the one-electron case, the description of the self-trapping 
phenomenon requires the inclusion of multi-phonon states. This is clearly displayed in 
Fig. 1^, where we have shown the weight of the m-phonon state in the ground state for 
various EP coupling strengths. One sees immediately that the maximum of \c^\'^ is rapidly 
shifted to larger values of m as Ep increases. Increasing the phonon frequency at fixed Ep, 
this tendency is reversed (see inset). In the extreme antiadiabatic limit {hu oo) we have 
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'rrimax = and the binding disappears for U > 2ep. 

As in the case of one electron it is interesting to look at the low-lying excitations of 
the inter-site bipolaron. Although we do not have a clear definition as to the momentum 
of this compound particle, it turns out that we indeed find a well-separated energy band 
if we again classify the lowest energy eigenvalues with respect to the allowed i^-states 
of our finite system (see Fig. P). The formation of the (inter-site) bipolaron band can 
be attributed to pronounced retardation effects [cf. the maxima in the nearest-neighbour 
correlation functions Cei^phi^) and Cei^eii'^) (Fig. 0) as well as the large bipolaronic hopping 
amplitude tb,e// (Fig. 6) at Sp = 3]. Surprisingly the dispersion of this 'quasiparticle' band 
becomes exactly like that of a free particle (with a strongly renormalized bandwidth) at Ep = 
U/2, where in the standard Lang-Firsov polaron theory the effective Coulomb interaction 
vanishes. As the EP coupling exceeds U/2, a deviation from the cosine-dispersion occurs 
and we expect that for Sp^ U/2 an extremely narrow AHBP-band will be formed. 



B. Holstein t J model 

Now let us turn to the case, where a few dopant-induced charge carriers (holes) cou- 
pled to lattice phonons move in an antiferromagnetic correlated spin background. In 2D, 



this situation, frequently described by the Holstein t-J model (^) ||5^j2^j25[| , is particu- 
larly interesting as it represents the basic electronic and phononic degrees of freedom in 
the Cu02 planes of the high-Tc cuprates. As yet, very little is known theoretically about 
the interplay between EP coupling and antiferromagnetic exchange interaction in such sys- 
tems. Of course, the exact diagonalization technique, as applied in the preceding section 
to the Holstein Hubbard model, provides reliable results for the ground-state properties of 
the Holstein t-J model as well. Here, however, one usually works near half-filling, i.e., the 
electronic basis is very large from the outset imposing severe restrictions on the dimension 
of the phononic Hilbert space. Therefore we are unable to reach the extreme strong EP 
coupling regime especially in the adiabatic limit. In the following numerical analysis of the 
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Holstein t-J model, the exchange interaction strength is fixed to J/t = 0.4 (which seems to 
be a reahstic value for the high-T^ systems). 

First, let us discuss the behaviour of the effective transfer amplitude, tp,e// = 



Ep^kini^p, J) / Ep^kini^p,^), shown in Fig. |T0|. Increasing the EP coupling at fixed phonon 
frequency hu = 0.8, the mobility of the hole is strongly reduced and an Holstein-type hole- 
polaron (AHP) is formed at about ~ 2.0. The continuous crossover from a nearly free 
hole-polaron (FP) to the AHP state is similar to that observed in the 2D single-electron Hol- 
stein model, i.e., at ~ a second maximum in the phonon distribution function (Ic™!^) 
evolves, which, for Sp ^ Ep, becomes more pronounced and is shifted to higher phonon states. 
For example, we get rrimax — 4 at = 4 and hu = 0.8. The increasing importance of multi- 
phonon states in obtaining the 'true' ground-state energy at large Ep becomes clearly visible 
in Fig. |10| by comparing the results for various phonon numbers M. There is, however, 
an important difference between the one-hole and one-electron cases which should not be 
underemphasized: In the single-hole Holstein t-J model antiferromagnetic spin correlations 
and EP interactions reinforce each other to the effect of lowering the threshold for pola- 
ronic 'self-localization'. This fact is in agreement with IMVLF-Lanczos results obtained 



recently by the authors As Fig. illustrates, the IMVLF-Lanczos technique, which 
variationally takes into account inhomogeneous frozen-in displacement-field configurations 
as well as dynamic polaron and squeezing phenomena, describes the qualitative features of 
the transition from FP to AHP states and gives a reliable estimate of the renormalization of 
the effective transfer matrix element tp^eff- Moreover, the IMVLF-Lanczos method yields 
an excellent variational upper bound for the true ground-state energy Eq, and therefore it 
provides an additional educated check for the minimal number of phonons one has to take 
into account within the Hilbert space truncation technique. 

By analogy to Eq. (^, we have calculated the corresponding hole-phonon density cor- 
relation function, Cho-phUi — j\) = (^ol^i^i^jl^o)' the 2D Holstein t-J model. Figure p!l| 
shows Cho-ph{\i — j\) as a function of the short-range EP interaction strength Ep at various 
phonon frequencies. The transition to the AHP state is signaled by a strong increase in 
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the on-site hole-phonon correlations which are about one order in magnitude larger than 
the nearest-neighbour ones. This indicates that the AHP quasiparticle comprising a 'quasi- 
localized' hole and the phonon cloud is mainly confined to a single lattice site. Increasing 
the phonon frequency the hole-phonon correlations are smeared out and the crossover to 
the small hole-polaron is shifted to larger values of the EP coupling. 



Now, let us consider the two-hole case. In Fig. |I2| we show the effective polaronic transfer 
amplitudes tp,e//(A^/i) vs EP coupling strength in the adiabatic {hu = 0.1), intermediate 
{huj = 0.8), and non-adiabatic {hu = 3.0) regimes. In each case we compare the one- and 
two-hole results to get a feel for hole-binding effects. Remarkably we find that tp,e//(2) is 
larger than tp e//(l) for ~ 1 and hu = 0.1, indicating a dynamical type of hole binding in the 
low-frequency weak-coupling regime where retardation effects become important. Indeed, 
the two-hole binding energy, defined as usual by E%{J,ep,huj) = Eq{2) + Eq{0) — 2Eq{1) 
with respect to the Heisenberg energy Eo{0), slightly decreases, i.e., hole binding is enhanced 
[E|(0.4, 0, 0) < 0], as the EP interaction increases at low EP coupling strengths. In contrast, 
at large phonon frequencies, with increasing Ep we find that E"^ increases, which seems to 



be an indication that retardation no longer plays a role ||2B|. On the other hand, in the 
adiabatic strong-coupling regime, where the two holes become 'self-trapped' on NN sites 
forming a nearly immobile hole-bipolaron, we expect an even stronger reduction of tp^effi'^) 
compared with tp ejj(l) (cf. the IMVLF-Lanczos results presented in Ref. [^. Here, a 
rather static type of hole binding is realized. 

To substantiate this interpretation we have calculated the hole-hole density correlation 
function 

CHo-hoi\^-j\) = mep,J)\h,h^mep,J)) (13) 

in the 2D Holstein t-J model. Note that Cho-ho{\i ~ i\) provides an even more reliable 
test for the occurrence of hole binding than the binding energy E\ ||5^. Indeed, when 
calculating E'^, we are comparing states with different quantum numbers, specifically with 



different S and S"^. In Fig. |I3| we present results for the non-equivalent hole-hole pair 
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correlation functions in the ground state of the Holstein t-J model with two holes. In 
the weak-coupling region the hole-density correlation function becomes maximum at the 
largest distance of the ten-site lattice, while in the intermediate EP coupling regime the 
preference is on NNN pairs. As expected, increasing further the EP interaction strength, 
the maximum in Cho-hoil'i ~ j\) is shifted to the shortest possible distance (remember that 
double occupancy is strictly forbidden), indicating hole-hole attraction. The behaviour of 
Cho-hoili — j\) is found to be qualitatively similar for higher (lower) phonon frequencies (see 
inset), except that the crossings of different hole-hole correlation functions occur at larger 
(smaller) values of Ep. In essence, our results clearly indicate that hole-bipolarons could be 
formed in the Holstein t-J model at large EP coupling. 

IV. CONCLUSIONS 

To summarize, in this paper we have studied the problem of (hole-) bi-/polaron forma- 
tion in the Holstein Hubbard/t-J model by means of direct Lanczos diagonalization using 
a truncation method of the phononic Hilbert space. Compared with previous treatments of 
the Holstein (Hubbard) model on very small clusters, we are able to analyze large enough 
systems in order to discuss polaron and bipolaron band formation, which has been a subject 
of recent controversy [0,0 . Our main results are the following. 

(i) In the case of a single electron coupled to Einstein phonons (Holstein model), we con- 
firm that the rather 'sharp' transition from a 'delocalized' nearly free polaron (FP) [or 
a large polaron (LP) in ID] to a 'quasi-localized' Holstein polaron (AHP) in the adi- 
abatic regime and the very smooth transition to a Lang-Firsov-type polaron (NLFP) 
in the non-adiabatic regime are both continuous. In agreement with recent exact re- 
sults [^,^,0 , we observe no ground-state level crossings or any non-analyticities as 
the EP coupling increases. We point out that in the one-dimensional weak-coupling 
case a large-size polaron is formed at any finite EP coupling. In the strong-coupling 
regime, the AHP state is characterized by pronounced on-site electron-phonon cor- 
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relations making the quasiparticle susceptible to 'self-trapping'. Most notably, the 
formation of an adiabatic Holstein polaron is accompanied by a shift of the maximum 
in the phonon distribution function to higher phonon states, which seems to be an in- 
trinsic feature of the 'self-trapping' transition. By contrast, the non-adiabatic NLFP 
ground state is basically a zero-phonon state. 

(ii) By calculating the spectral properties of a single electron, we have found convinc- 
ing evidence for the formation of a well separated narrow polaron band in both the 
adiabatic and non-adiabatic strong-coupling regimes. In addition to the expected 
band-narrowing we also found a deviation from the 'cosine'-dispersion away from the 
adiabatic and antiadiabatic limits. Although the 'coherent' bandwidth, deduced from 
our finite-lattice ED data, becomes extremely small in the adiabatic strong-coupling 
case (polaronic band collapse), we believe that the AHP does not lose its phase coher- 
ence and can move itinerantly. 

(iii) Investigating the two-particle problem in terms of the ID Holstein model, we could 
clearly identify the transition from a extended (large) bipolaron to a 'quasi-localized' 
(on-site) bipolaron (AHBP) as the EP interaction strength increases. Stabilizing a 
two-polaron state in the weak EP coupling regime by taking into account the on-site 
Coulomb repulsion (Holstein Hubbard model), we found a transition to an inter-site 
bipolaron at about Sp ~ U/2. It is worth emphasizing that this inter-site bipolaron 
appears to have a dispersion that resembles very closely the cosine-dispersion of a 
noninteracting particle with a renormalized bandwidth. If the EP coupling is further 
enhanced (^p ^ f//2, Tiuj, t), a second transition to a 'self-trapped' on-site AHBP will 
occur pT| . 

(iv) Analyzing the hole-polaron formation in the framework of the 2D Holstein t-J model, 
we found that the critical EP coupling for the polaron transition is substantially re- 
duced due to 'prelocalization' of the doped charge carriers in the antiferromagnetic spin 
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background. Therefore we suggest that polaronic effects are of special importance in 
(low-dimensional) strongly correlated narrow-band systems like the nickelates and 
high-Tc cuprates. 

(v) Regarding ground-state properties of the Holstein t-J model in the two-hole sector, a 
detailed study of the hole-hole correlation functions and the two-hole binding energy 
was carried out, yielding strong evidence for an enhanced hole attraction and the 
formation of hole-bipolarons as a dynamical effect of the EP interaction. 

Of course, the exact results presented in this paper hold for the Holstein Hubbard (Hol- 
stein t-J) model with one and two electrons (holes) on finite ID (2D) systems, i.e., we are not 
prepared to prove any rigorous statements about the thermodynamic limit here. However, 
we believe that our main conclusions (i)-(v), in particular the existence of well-separated 
polaronic and bipolaronic quasiparticle bands even in the adiabatic strong-coupling regime, 
will survive in the infinite system. 
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FIGURES 

FIG. 1. Ground-state energy Eq and weight of the m-phonon states |c"*p (inset) as a function 
of the maximal number of phonons M for the ID single-electron Holstein model on a four-site 
lattice. The model parameters are: £p = 6.0, hu = 0.4 (all energies are measured in units of t). 

FIG. 2. Effective hopping amplitude, ip,e// vs Sp, for a single electron on a ten-site lattice 
described by the Holstein model. 

FIG. 3. Electron-phonon density correlation function Cei-phHi—j]) as a function of the distance 

i — j and EP coupling £p (inset) at huj = 0.4. The results are given for the single-electron Holstein 
model on a eight-site chain. 

FIG. 4. Single-electron band dispersion Ek — Eq at £p = 3 in units of the coherent band- 
width AE. The inset shows the distribution of the eigenvalues at hu = 0.4, where the horizontal 
axis counts the eigenvalues sorted by magnitude and the vertical axis gives their absolute values 
En := - Nhu/2. All results are given for the eight-site chain with M = 18 phonons. 

FIG. 5. Band dispersion E{K) of a single electron described by the Holstein model on ID rings 
with N sites, where Sp = 3.0 and huj = 0.4. The insets show the low-energy part of the one-particle 
spectral function Ak(E) taken at the if -values indicated by arrows. The dotted line corresponds 
to the dispersion of a free particle with a renormalized bandwidth. 

FIG. 6. The effective polaron and bipolaron transfer amplitudes, tp^gff and tb^eff, are shown 
as a function of the EP coupling strength Sp in (a) and (b), respectively. The results are given for 
two electrons on the eight-site chain with M = 21 phonons using periodic boundary conditions. 

FIG. 7. Electron-electron [electron-phonon] correlation function Cei-ei [Cei-ph (inset)] vs Sp 
at = 0.4 for = (a) and U = 6 (b). 

FIG. 8. Weight of the m-phonon state in the two-electron ground state of the Holstein Hubbard 
model, where U = 6 and huj = 0.4 [hco = 3 (inset)]. 
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FIG. 9. The lowest eigenvalues Ek for the two-electron case on the eight-site chain, where 
huj = 0.4 and U = 6. The energy band is given in units of the coherent bandwidth AE = 0.0481, 
0.0217, and 0.0190 for Sp = 2.5, 2.75, and 3.0, respectively. 

FIG. 10. Effective hopping amplitude tp,e// and ground-state energy Eq (inset) as a function 
of EP coupling strength Sp for the 2D single-hole Holstein t-J model at J = 0.4 and huj = 0.8. 
Exact results for the ten-site square with different numbers of phonons M are compared with 



approximative IMVLF-Lanczos data |24]. For further explanation see text 



FIG. 11. Dependence of the on-site and nearest-neighbour (inset) hole-phonon correlation 
function Cho-phUi — j\) on the EP coupling Sp and phonon frequency Huj at J = 0.4. Results are 
given for a single hole on a ten-site lattice with at most M = 12 phonons. 

FIG. 12. Effective transfer amplitudes, e//(-^/i) £p, are compared for the one- and two-hole 
cases at various phonon frequencies, where J = 0.4 and M = 12. 

FIG. 13. Hole-hole density correlation function Cho-hoil'i ~ j\) of the 2D Holstein t-J model 
shown as a function of the EP coupling Sp for J = 0.4 at huj = 0.8 and huj = 0.1 (inset), where the 
numbers 1-3 label the non-equivalent distances — Rj\/a = 1, \/^, and \/5 between NN, next 
NN, and third NN sites on a ten-site lattice with M = 12 phonons, respectively. 
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